nil -a 565 5 


NASA CR-14180 



SPANWISE VARIATION OF 
POTENTIAL FORM DRAG 


W. C. CLEVER 

Los Angeles Aircraft Division 
Rockwell International Corporation 
Los Angeles California 

MAY 18, 1977 


Prepared for 


IVJ/NSA 

National Aeronautics and 
Space Administration 

LANGLEY RESEARCH CENTER 
HAMPTON, VIRGINIA 23665 


MAY 1977 

RECEIVED 

NASA STl FACILITY 
INPUT branch 

b'nc# 





SPANWISE VARIATION OF 


POTENTIAL FORM DRAG 
By W. C. Clever 

Los Angeles Division, Rockwell International 


SUMMARY 

A new method has teen developed for calculating the spanwise variation 
of potential form drag due to thickness and lift for a nonplanar wing of 
arbitrary planform in subsonic and supersonic potential flow. A computer 
program has been developed to perform the numerical calculations. 

The configuration is subdivided into a large number of panels, each- 
of which contains an aerodynamic singularity distribution. Linear ly -varying - 
chordwise source and spanwise vortex distributions are used to represent 
thickness and lift. The normal components of the velocity induced at 
specified control points by each singularity distribution are calculated 
and make up the coefficients of a system of linear equations relating the 
strengths of the singularities to the magnitude of the normal velocities. 

The singularity strengths which satisfy the boundary conditions of 
tangential flow at the control points are determined by solving this system 
of equations. 

The potential form drag due to thickness and lift is determined from 
the strengths of the singularities and the panel drag influence equations 
obtained by analytically integrating the tangential and normal velocity 
perturbations over , the panel" surface respectively. 

Several examples of the spanwise variation of drag calculated by the 
program are presented and compared with available analytic results . Good 
correlation was obtained in each case. 
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INTRODUCTION 


Prior to the present study, no procedure existed for computing the 
spanwise variation of potential form drag at transonic and supersonic speeds. 

Also, at subsonic speeds the calculation of the spanwise variation of drag 
was limited to planar wings represented by a vortex lattice with even spacing 
in the spanwise direction. Due to the necessity of determing the leading edge 
force in order to obtain the nonlinear vortex lift and drag on nonstraight 
leading edge wings (using the Polhamus leading edge analogy), a solution 
to this problem is essential. Also, in order to predict the contribution of 
horizontal surfaces to such stability derivatives as , C n an< ^ C n 

and the pivot moments due to the outer panel of a variable sweep wing, the 
spanwise distribution of leading edge force must be computed-'; 

In addition to the requirements for the spanwise distribution of leading 
edge force for the prediction of loads, the analysis is also valuable in optimizing 
the planform shape to minimize three dimensional induced effects. For 
instance, the occurrence of local wave drag, shocks, vortex formation, and 
flow separation induced by the planform shape can be minimized by shaping 
the planform to have a constant value of leading edge force across the span. 

In order to compute the section potential form drag for wings of 
arbitrary shape, a finite element must be used which is capable of representing 
the waves produced by supersonic leading and trailing edges and which has a 
spanwi se variation of vorticity that is continuous and goes to zero at the tips of 
the element. The most elemetrary lifting finite element which meets these requirements 
is one in which the bound vorticity varies linearly in the spanwise direction, 
such that the vorticity at the tips of the element is zero, and is constant 
in the chordwise direction. The influence equations for this type of panel are 
in terms of the same quantities as needed in the calculation of the influence 
of a constant pressure panel. Therefore, the linearly varying pressure panel 
should cost the same as the constant pressure panel to use. . However, due to the 
lower order singularity at the tips, the normal velocity can be analytically 
integrated over the complete element to obtain a finite value for the potential 
form drag due to itself or another element, a finite element with a spanwise 
constant and chordwise linearly varying source distribution was chosen to 
represent the contribution due to wing thickness. The chordwise linear 
variation was necessary to avoid the infinite drag associated with constant 
source panels having sonic leading edges. 

This report describes the mathematical development of a method for 
calculating the spanwise variation of potential form drag of a wing at subsonic 
and supersonic speeds using these linearly varying panels. The wing may be 
of arbitrary planform and nonplanar provided the wing panels are all parallel 
to the aircraft axis. 
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LIST OF SYMBOLS 


b 

B 


c avg 

C D 

°Dij 


CP 

ACp 

Cr(y) 

k 


R 

Sref 

T 

u,v,w) 

Vr,V e 

P > 2 

* 

7 

$ 

§ 


Span 

Chord 

Reference chord such that b x c a vg = S re f 

Drag coefficient, Drag/q „,S re f 

Near field drag coefficient giving the drag 
singularity j induces on area i for unit values 
of pressure coefficient 

Pressure coefficient (P-P* )/yi° 0O 0^ 

Net lift pressure coefficient CPiower-CPupper 

Leading edge thrust coefficient, Thrust/unit length 

Basic functions integrated in Appendix C 
Vector of basic functions f * (See Appendix D) 
Coefficient = 1 if 1. or 2 if M „ > 1. 

Free stream Mach number 

X l 

Free stream dynamic pressure r ~ 

Either V x% or sj , as indicated 

Wing planform area 

The tangent of the sweep angle 

x, y and z components of velocity 

Radial and circumferential velocity components in 
polar coordinates 

Either -1 or 1 - 

Usually X-T^ 

Usually T 

Usually Te 

Perturbation velocity potential 


THEORY 


FORM DR/VG DUE TO LIFT 
Gene ra 1 


It will be shown that an <'xact solution for calculating the 
form drag due to lift can be obtained by integrating the product of the - 
local induced normal velocity times the local ACp over the planform area. 

If this integration is performed analytically, the drag 
calculation for that pressure distribution is obtained exactly, and the 
drag as calculated in the near and far fields must agree. This drag 
includes the leading edge suction forces since all momentum losses have 
been accounted for in the far field analysis. To be suitable for analytic 
integration, the spanwise distribution of vorticity must be continuous over 
the planform, so as to avoid trailing vorticies of finite strength and the 
resulting infinite drag. Any method which does not employ a continuous span 
wise distribution of vorticity can only calculate the drag approximately 
and the value as calculated in the near and far field need not agree. 


The method employed in this effort uses the most elementary method 
of obtaining a continuous spanwise distribution of vorticity. The planform 
is composed of quadrilateral panels each having a linearly varying vorticity 
distribution in the spanwise direction and a constant vorticity distribution 
in the chordwise direction. Two adjacent panels are combined into a finite 
element having a vorticity distribution which is zero at the outer edges of 
the combined panels and varies linearly in between, reaching a maximum value 
at the junction of the panels. 



The entire vorticity distribution can be described by knowing the value of 
h Cp at the junction. A finite element extends over two span segments 
of the wing. Therefore each -span segment has a distribution of a Cp which 
is the sum of the <dCp from the ..left, panel of one finite element and from 
the right half of another finite element. On span segment L 


AC p 


AC e. 

«.♦ I 




+ AC p 

l 
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The influence equations for the type of panels required for this distribution 
are developed in the next section. The values of ' Cp for each • singularity are 
obtained by satisfying the boundary conditions for normal velocity at a set 
of control points distributed over the planform. The control points are 
located at the same places used for constant vorticity panels. 


The section after that describes how .the normal 
velocity induced by one of the panels making up such a panel pair is in- 
tegrated over a panel of another panel pair after multiplying 
by the local a C p. Since the integration is performed exactly the result 
is the drag which one panel induces on the other. Since the drag is 
proportional to the product of the aCp's at the panel edges, a drag 
influence coefficient Co^ may be defined such that the drag which panel j 
induces on panel i is 

Cp ■ • ACp. AC P . 

t 4 



AC.p 



induces a velocity normal to 



segment i 



i 


The integral of induced downwash times A Cp is - 


Cp,. ACp.ACp. 
‘4 * ■* 
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Since a finite element is composed of two adjacent panels, the drag 
which a finite element (composed of two panels) induces on a panel located 
on a given span segment can also be written as * 

Co ACp £iC», 

i.1 ‘ •* 





Since a span segment is composed of a iCp from two singularities, 
the drag induced by singularity j on span segment i can be witten 


ACp, +C, 


AC, 


I*',, 




AC P 

j 


The superscripts L and R refer to the left or right half of singularity, 
pressure distributions on segment i» The drag induced on a panel spanning 
such a span segment is obtained by summing over all panels j . 


The leading edge suction force is obtained by subtracting the near field 
drag, as described above, from the zero suction drag. On a given panel the 
zero suction drag is calculated by integrating the product of the local ACp 
and the normal velocity as given by the camber shape. This normal velocity 
js jonly equal to the induced velocity at the control point and the difference 
accounts for the leading edge suction. 
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Spanvise Linearly Varying Vorticity Finite Element 


A constant pressure or constant vorticity panel with a quadrilateral 
shape can be constructed by adding or subtracting four semi-infinite tri- 
angular shaped panels. These semi-infinite triangles, each determined by 
a corner of the quadrilateral, can be assumed to induce a velocity potential 
everywhere in the flow. However, each corner represents only an integration 
limit, and all four corners must be included to make any sense. If it is 
kept in mind that four corners must be included, one of these triangles, 
having unit ac p , lying in the E=0 plane, and having sides determined by 3=0 
and X-Tj= 0 , induces the following velocity potential 




8tt 




Ta r 

T 3 (x-T a )-T*i‘ 




where 


- 2/TV7JL0 


/rV ft + (Tn*j a*a) 




1 


/= \-n] m 



/&* >0 
/&*<o 


To obtain a panel with a spanwise variation of AC P we can combine 
panels composed of thin chordwise strips, each with its own ac>, and then 
take the limit as the strip width goes to zero. 



The resulting potential is 




- J -^rtooL*- Ll s> * 

• >* 
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s'uch 0 ^ 86 Ve WlU ** a l inearIy 


trying Cp( Q an(J 


Therefore 


c * ( 5>= o,>u- 

* b a US) - T ? c .,-, 

3 c p(a) = b 


straight leading 


*<*>3, 2) - o. 


3^ 2) 


vhere^if We let 


(a ' +b Wt,[(x-T ;i ) 4T(3 } 

'»*• 3 J 2- J 

tb } r 

* w 


*»ak 


j = 


X"T» 


T <a-3) 

I * Ta 


8 = 




the " tI<*-T a >.T( vS) , <3 . 3)( J = 


*£<-t s >.t (v3>j a .j >S J 


I)ct 


7 


r ^ 

Tj ^(S»7,J.7,^ 

T<a -=W 

t(i,7,1,T)j 3 

C£ <x ' T a), T <»-5> >T * >T j| a 

$• ^ 

*4* 

*e/ x > 3, i) - * r,* _ . 

13 w, v^, ij 
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However >the integration over can be performed more easily by 
integrating the velocities rather than the velocity potential. In our case 


*J*>a>*) = 


A 


8ir 


> i o R+t***?) * • o*, BR lUt aN) 


( t ♦ a'-?; 




and 


where 


U oo S 

3^o 

3X 

3<>o 

35 

A 

8TT 


v 0e * 


■I 

3F 

aj» 

- ‘ 

* #p 






<0 = 

34. 

32 

- T -£*2. -- 

3 5 

At 

8ir 

J i 8._ (l+lm) . 7R 

^S S * 


2 ^JBR- (5te*l) 




R' - (5**7)*+ ♦**) 
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If we define 


v ol (u^&) 
B) 


j u oo Cl,7,^,6>«»7 

J V oo Ci,7,l,B)J7 

j w oo (s i 7,^0)«*7 


then 


u.„(j,i,s,b> - -|r{ * l ir-h, 


QR.+ ii-a'V 
a r. - 1 5 + Q* 7 ) 






At 

0TT 


j (7~S) *<>~- 


-15- + 

75-i l 


iC8V.^i 


aRt + 
6R.-C5 + a'-7) 


Txw! - ? - z 5 i V‘ s,) ]} 


Ar f , ft R.+ (5 + ‘7> /v * . ' • 0 GR + (Sf0j7) 


+ l**-” ^T* + 2 * 


- (2-A>[(7-5>Yl<^C7nV 21^* Y - 2 (5*7>] ■ 


These velocities can be thoucht of as being induced by a semi -infinite 
vorticity panel with C P t^) - =• 7 
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For triangular loading, consider two adjacent panels in the z=0 plane. 





Then for any velocity component, say w, we can combine terms to obtain 




*) = Cp | - W^x-x.jVi, J 2, 


T o) + wj X-X, , > T, ) 


+ T -*,(***■> 3**. >*> T o>] 

o 


- -7- :[ U V< *-*•» V*.» *> T . > - *>T ( )] > 

T, '•Oi' V J 


T o <W 




T 0V»,> 




j UJ^ X -X 'S * > 4 , , ^ 


where u) 0o is the velocity induced by a constant vorticity panel and u) ol 
is the velocity induced by a spanwise linearly varying vorticity panel. 
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Interference Drag Between Two Spanwise Linearly Varying Vcrticity Panels 


Consider the following definitions 

S = *-t^ 

1 

S ^ r a 

R* = C? + V) 1i+ Co 4 -OC7%i*) 

Then, .considering only the near field terms, the velocities induced by each 
corner of a constant vorticity panel are : 


IX (SUB) s 

DO > ■*) J 


4 -| SR. 

/AA. _ 

Sir S7-S* 







T+x*J 




8 tt 1 2 ^ »va» 7 ) © 2 aR-dta’ 7 ) + 7 * 4 i l 


A. 


t / 5 l > O , n_<| 

z y < o , "»->> 


Wow define for each velocity component 


JJ u i; U,l,I.6) 

d 

U>. . / 5 7 ^ qs 

^7 c 4 v ' '» ' J 




A discussion of the calculation of these velocity components is 
given in Appendix C 
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Using this notation the velocities induced by a single corner of a 
spanwise linearly varying panel would be : 




because 


u. * u. 

a -» •• 00 


— v » v 

0| 00 


37 


cO * u? 

9*7 °» oe 


The factor — is necessary because the integration is over 7 and not y. 

r 

Using the above definitions the interference drag between two panels 
can be derived. Consider two panels A and A having leading edge sweeps T and 
T, and coordinate systems such that A lies in' the plane Z = constant and A in 
the plane Z = constant, and the coordinate axes x and x coincide. 


A 



Ca-V 

U-V 
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The coordinate transformation between the two coordinate systems can be 
written in the following form 1 

A A 

X = X X - X 


a -- 


* 

-/(.'i-tS 




3 


i 

/(l *Z l ) 


IV-ri) 


z 


I 


(- h) 


z = 




(*V z) 


If we define 


T -/l + X V 

CA + i) - — 

A 

t -- - r 


(JUi) 


T -/l+T* 

A 

T 


then the following relations, which will later prove to be useful, can be 
derived : 


(JUOCA+i) 


0 *t x ) 


Co x -Q _ Cs x -Q 

U*i) " (i+O 


A - - 

a 

<A*l) 

(A-t x ) = 


(A-f‘) 

Ca + O 


A * t 

U-T ) 
cii-f 1 ; 

l(m‘) 
Ca* 1 ) 


(A+i) 

Cast's') 


(.A ♦ 1) 
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In addition 


S = * - T £ 


7 




(3»”) 


A 

l 


A A 

T Z. 





u»o 


c 7 + r j ) 


U»o 

0+T\> 


U- 77 ) 


and 


( £ * 7 ) - Cl + 7) - £ * x 

UUiHtM*) - OUOt7**j‘) 

R* * (5*7)* •*• CQ l -0t7 + S*3 


(.1*7) * 


tB v - 0(7** 


I 1 ) 


Assume each panel has a vorticity distribution 


ac p= ^( 3 - 3 .) 0,0 A 

^ c e • *,C3-V - A 


Then the velocities induced by the leading edge of panel A at any point 
(i,7,$) are: 


u. 


V 


to 



C 1-1,, 7-7. 


5-5.,S) - v .,^-l t ,7-7 tJ i-S 0 B^- (7 r 7 ( )V e (M t ,7-7 4> S-J„B)J 
- v0 Ci CS-S,,7-7 li S-S, J B)-(7 4 -7,)*O o .U-I J ,7-1 ij S-S lJ B)} 


Note that 5,= i jt because (.x^-t^^cv 7 !,) on the leading edge , and * 5 t . 
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The induced drag of panel A on panel A is obtained by integrating the normal 
velocity ( 1 * 1 *)'^ (u>-rv) times the vorticity strength ac* » 
over panel A. feee page 34.) • At a fixed value of 7 the integral in the x 
direction may be peformed at once. Since uO a to we can write for 

the integrals Of and u»: ** ** °' 


C [ l' 

*T.t. 

*r.«. 


* and 


where ^ L.E.^and ^ T.E. refer to the values of 1 at the leading and 
trailing edge of A for some value of 'f 


Since the procedure is identical for determining the location of both leading 
nnd trailing edges, only the leading edge - will- be considered, -Mone this edge: 


<*- > 


* t •/• ♦ t * (a- * 3 :> 


U-e.) - <**-**:> - * cz 4 - z t ) 

S- x [Vui 1 t-t]( v3i ) 

^ CV + Z L 1 't»> 
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To complete the integration over a quadrilateral panel we insert the 1 limits 
(which depend upon ^ ) and integrate over ^ . 


If we define 




the induced drag is then obtained from the following expression. 


Cd S 

c 


% 


* — J O-V ^ U V V * «1-V> 1-1, . (i,- 1,) » 8]^ 


t 


O+tM . ^ , 

- 73 J <V s 4 > tr <7*V>°] J 7 


7, 


c. + ts r 


7.1 


These integrals require the following quantities of type W such that: 


v/ - 7 3 .. 

v,’ - 

where 

v ' * Jw - y-o-^-v 
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Integrals of this type are discussed in Appendix D 


Since 


i - 5 Cf ) and \ - SO?) , 


we will denote 


■k by 4..C7) 

‘A 1 J ‘J 


Therefore, since 


-.'hero 


fc* 


A = 


sir /iTP 


A 4 -..ct 4 


-5 

Z 

Z 


we have 


8TT Co S « t , 

* C.c p 


7 

j< 

* •A 




- j 


S,&>] 

l /&>o 
i /2>'<D 



Or, using the results of Appendix D, 


8rr Co -Skip 
* C P C P 


= ^T~ a |(V7>ElVvi)- F L l V^]- E *t, ( V^ + L EP * FE ] £ ?JV’0 

- (V 1 )£ [ V £*f ( 1 3 - 1 ) - [ EP - PEje fjvi) 

l 7 ’ 7, 

- n- 7 a )El? ti ci- 7 )-F? je tv^ 

> 7 , 

- AE ( 7 t - " 7 ,) | - Eljvv 


- <Wf t /W + E UW + t V^>lOv 7 l )-T % cvv]J 



U*l) 

0*0 


A 


{ 


(VI) plT ^ 0 C V 7 >- 


U*»l 

0 *t l ) 


F V + A ,( V 7 ) - l P E ♦ E F ] F T T^ o ( V 1 ) 


-(. 1 ,-DF'r ^V 7 )*^/H,<V^ + l F6 * £F J FT L ( V^ 

, 7 4 , ,> 1 , 

♦-C 1 -V f1t l i ( V' 5 ) ‘ VV 1 )][ 


^vvAj-tw^OVv* 


+ a-\) FT L ( V.v - £FT l ( vi)-f^ ) FT Ovv 



sir 


Co -Sate 


Jk 



'/\+T x Ovt*) 

t 1 tl*D* 


Now define 

AS 

tt 

and note 


c P 2 . p 

(E P ♦ PE) E ?J1 4 - V , ) - e‘ T a C 1 4 - V + <V 1 A ) E [ ft V 1, > - F f ao C V V , )] 
( e f + pe) e 


■ (eF * FE > £ OvV-^ Sl tVV - I 

■ uf+ Fe)E ^ c v^> +E ^ 1 t v i « ) ~ ( vv e tvvv- F L t v^>] 

AT-| (EF*FE)FT? 4 C 0 \- 7 ^^/^ 

- (E F ♦ F E) FT ^ 1.) - ^ F 'T ? 41 C V 1, ) 


CE^pejFT^-v^^pV^vi.) - tVV F * T L< V V 


i EF * FE) FT FT UVVHI-^PjFT TjV V» EFT? <1,-^] 

+ tv V F l FT t*t vv - ( VVT +VVV] } 


<W ' T l W 


AS 




TFT 

t"f 1 t 


- CA*0 E 


A. 1 A 1 

a*o e 


T E T 


C&*|) 

(J*r*> 


T EFT 


isUo* £ * 

0 *T*) he 


t''(ef + F£)FT v CEF-F£>e 

ci+t*; 


19 



Therefore 


Bn C p 

A 




/TrT 1 


r~ A J==r- ElVVV*f?A^] 




+ («*«)£ f 4> ci- 7.) -e‘T ai (7i7 t ) - 


*T*AS — 

(« . f«i t • £‘ 4 „(VV • Jj# t[ Wv- '<JVV) 

■>' e I E Ow ■ yiTft 1 + ,. l -vv] | 


+ T £>& 


| The same terms with the following replacements 


A 

1 — 5 

A 

7—7 

A 

5—1 

A 

e> — e> 

T — o. T 

A 

A — A 

A 

r —* z 


I — 3 

Z — A 

3 — i 

A- — Z 


Prom this expression it can he seen that the total interference 
drag between two panels is zero. Interchanging the influencing and 
influenced panels changes the sign of the expression and the sum is 
zero. 
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There are four expressions (AE, AEE), AEF, and A(EF + FE)f required for the 
evaluation of the interference drag. These coefficients will now be 
evaluated . 


AE 


= AH 


CANtV) 


~\ 

a-z x ) - jr l (s l : i) 


- a'u-O 


a-r 1 ) >■ 

T (a' a) 

za*<) 


E -- H 
F - H ! Ci 


AF . A h'g, 


(A* 4 T l Q*) 


— . * 




U-T l ) ~ 

- Q 1 C A - r* - ) f r'tB l +A) + 

a-T l ) - r l (A + .) 
t[(B%A) +UL-T 1 ) + ^T*Ca l - 0] 
T [ (£+•) * (A-T 1 ) ] 


>- 


AEE = f\ h" H *' 




JUA-t*) - - ir\8-0 


- B l [ Ma-t 1 ) - t'cb’va) ] 


■A T l ) - t \ &%*■) 


T ( A CB^A) + Q l (Jl-r 1 )] 
T [ AU-I) '♦ 
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A F £ * A H'g, H 


I 


AEF 






- + t*(&+«.)*V jr\iV rV)(a-o 


(JL-r')*- t‘( u**) 1 * Cb >0 ♦ ***«*•)] >- 


z*LSL-t x ) ($.* o - j r*(a-o 




r 


ah" HQ 


^ * 


(A - 1 * ) - 1 * [(»♦•>*■» { I ♦ <*♦») - i OT *) }<a*- 1 ) ] 




- t 1 -*[i- 

(JUr 4 )*- r‘ [<*♦»>*♦ cb'-o ^ 

2TC4-t*)CoV«.) - jU*0) 


2 rU-t‘)lA *0 - yT*ta-o 

3U 


A[F6^EF]e • A [h'<SH''« H’V<0 H"’ 

(A-^) 3 +r‘t(i*r*j-w*<*- ~ »5iouo - 2 t\b*-o} ( o'-oj 

2 , 

- AU-t\)V ♦ T l JUQ*J.) r tzVt*-T')(tt*A)-± [f -<4 -m)] U'tV^b’-O 

— <! U-t l V*+ x* [o+t‘>*-(a* 0 4 - { * ♦ aCA-T*)- jU«.o] 1 Cs^-o] > 

(A\tVy* I L J J ( 

r| CA-T')V- T\q\s.^* 2JU*-r*)(B + 4) - ^ C4*»T‘&*} r l (a l -0 j 
r{ U-r 1 /- t x U*‘)'+ eA(«-r*XA*o •»£[ i- £(*♦«)♦ j t »♦**>] j 
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In terms of more familiar quantities 

£ z — 1 — (t - t ) * b 1- * ~h: (t 1 * a )' 





(jC+t'q?) 




|(t - T c&4<f> ) + A***- ^ i T * A ) ^ 


. Z Z , 
C©-^ tf> 


| (T-r> 1 


z (l-C^f)TT 


b* ^ 


(^5.- t l ) - -1 (.T - T CLO+-4 ) Cer^d) - T \ - ( T C - T ) 

T cWf L j T <K^</> 


£jO , — — ! 1 ( T V/i^ CoA<J> ♦ TCT-Tc«^J » | TT + A* | 

T Cs&c±*i> ' ^ 


I *01 

T Cp 

T 


OO 


T coc> tf> 



23 



^AE*- 


3 


T(T-Tt^f)(T^-T)- J^y[T(T t + / &)+j-T/b Z C**t] 

- — (T^/S*) [ (T-Tc^^iTCt^tf-T)- JusSj (TT + /J> l c^A^ )] 

T [ (T-ttM^)(f CJB\<t>-T)- (TT >- 

[(T- T c*a <f> )( TT + /J 1 tM^)* (T c++<f -T >(T +y& 1 ) ] 


Ta !~4 [ (T-T)(t*t)”y / 2>t ■***} f ] 
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For the case where JL, = 0 or T = T and 2 is constant, then. 
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The expression for Cj) with t.t c^n be shown to be equivalent to this 
in the limit as t-*-t. This can be demonstrated by using the fact that 
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FORM DRAG' DUE TO THICKNESS 
Chordwise Linearly Varying Source Finite Element 


influence functions for a constant source panel are 2 
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Integrating in the chordwise direction gives the influence equations 
for a chordwise linearly varying source panel * 
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The source panel singularities have the following properties: 

1. All panel leading and trailing edges are at constant (^-); side 
edges are at constant ^ 

2. Each source finite element is composed of a panel pair 

3- Source strength varies linearly with chord measured from the 
leading edge of a panel pair, i.e. the maximum value of the 
source strength is proportional to the local chord and attains 
this maximum on the panel edge joining the panel pair 
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The panel pair induces the following perturbation velocities 
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Results using these equations are compared with 
an exact solution for a biconvex airfoil in figure L. The only- 
difficulties encountered were due to the fact that a chordwise linearly- 
varying source panel must have zero thickness slope at the leading and trailing 
edges. Therefore the thickness variation could only be satisfied approximate' •• 
in these regions. 


Interference Drag Between Two Chordwise Linearly Varying Source Panels 


To obtain the interference drag, the u. velocity induced by one finite 
element must be multiplied by the local slope on another finite element and 
integrated over the area. 
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Figure 1. Rectangular Biconvex Wing AR = 5 Pressures at Mid Span 







Consider first the integral over i of a typical term 
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However, analogous to the vorticity case, 
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Performing the integration over 7 of a typical term (where S-Kt) , x> - Ui) ] 
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Therefore, one typical term for the drag induced by one panel on emother is: 
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- the same term with the indie ies 1 and 7 interchanged 
and all variables based in the (M,3) coordinate system. 


The complete drag coefficient requires performing this calculation for all panel 
corner pairs and adding the result- after mutliplying by the appropriate sign. 


From this expression it can be seen that the total interference drag 
between two panels is zero. Interchanging the influencing and influenced 
panels changes the sign of the expression and the Siam is zero. 
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In the foregoing 
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NEAR FIELD AND FAR FIELD INTERFERENCE DRAG EQUIVALENCE 


This section will show the equivalence between the drag force 
as calculated in the near field and as calculated in the far field. The near 
field drag induced between two singularities is calculated by multiplying the 
induced velocity normal to one vortex panel by the local Cp and integrating 
over the area. The far field drag is calculated in two parts by integrating 
perturbation velocities over the appropriate areas on surfaces far removed 
from the singularities. 

The singularites chosen to show the equivalence are two infinitesimal 
areas of constant pressure panels. It will be shown that the total 
interference drag between these two singularities is the same whether 
calculated in the near field or far field. The velocities induced by these 
singularities may be found by differentiating the velocities induced by a 
constant pressure panel with respect to x and Y. These singularities actually 
represent elementary horseshoe vorticies. Any singularity field or pressure 
distribution on any set of surfaces may be constructed by integrating 
together the appropriate distribution of these elementary elements. 

Consider two infinitesimal area vorticity distributions and their 
associated trailing vorticity. Place an (x, v, z) or (x, r, 0) coordinate 
system inside the upstream area with the y axis passing through the trailing 
vorticity. Then later in this section, the following results are derived 
for supersonic flow. 
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The near field drag is 
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the far field wave drag is 
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the far field vortex drag is 
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P»*on the''-!' results it is easily um n that the total far field drag is 
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For- supersonic flow only the near field terms of the influence equations 
were considered.For subsonic flow the near field drag as calculated using 
these terms will cancel since the value of x will change sign when we 
interchange the influenced and influencing vorticity distributions. For 
any distribution of vorticity this interference drag will cancel exactly if 
the integration is performed exactly. Therefore in subsonic flow the only 
contribution to the total drag comes from the far field terms of the 
influence equations and this is identical to the vortex drag. 


Induced Velocities - Supersonic 

To derive the above equations consider a constant vorticity panel 
located in the Z = 0 plane with one of its corners at x = 0, y = 0. 
Consider also a polar coordinate system with the following definitions 
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The velocities induced by an infinitesimal streamwise- running strip of this 
constant vorticity panel are found by differentiating the velocities induced by the 
edge of the panel at y = 0 with respect to y. The resulting strip is actually a 
line doublet which varies linearly until the end of the vortex strip and is. constant 
thereafter. The velocities induced by the end of the strip starting at x = 0 are: 


u. = 

3 

IA = - 

9^ 00 

c p 

4-TT 

X 2 : 
v'R 


Cp*a 

4tT 

- 

v -- 

_ J_ y _ 

3 & 00 

Cp $a 

4-rr 


[ y&V 1 ] 

4tt 

t«VVJ 

u) = 

d 

- — lO - 

oo 

S 3 
4-t r 

{■& 





and radial and tangential velocities in polar coordinates may be defined 
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The perturbation velocities are found by evaluating the functions at the start and 
end of the strip using a coordinate system located at the respective end points. 


The velocities induced by an elementary horseshoe vortex or line doublet are found by 
differentiating these velocities with respect to x and multiplying by -1. 



V 


Cp&S 

^ ° f 

4ir 

4tt R!* 

Cp&S 

X 

4TT 

r l R> 

Cp&S 

x 

4TT 

r*R* 


* &x 


36 



Near Field Drag 


Consider an elementary horseshoe vortex with its local z axis inclined at an angle (f> Q 
vith respect to the vertical. Then the velocities induced are 



If we consider two singularities, denoted by subscripts 1 and 3 then the 
component of velocity of vortex 3 normal to horseshoe vortex 1 is * 
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Since only one singularity can lie in the downstream Mach cone of the other 
the near field drag is given by this expression on page 35* 
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Far Field Drag 


The far field drag is claculated by integrating over two surfaces at infinity? 
The wave drag involves an integral over a cylinder whose axis is the x axis and whose 
radius approaches infinity. The vortex drag is found by integrating over the end of 
the cylinder as x -» <» (Trefftz plane). 
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Wave drag -to evaluate the wave drag we must perform an integration over 
a large cylinder whose radius approaches infinity. To avoid complications 
in the integrations we will calculate the interference drag between two 
strearawise vortex strips of finite length. The wave drag between the 
elementary horseshoe vorticies may then be obtained by differentrating the 
result. The velocities induced by these finite length vortex strips are 
given by the equations on the top of page 36. 
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Consider a coordinate system located at the center of the large cylinder. One end 
of one of the vortex strips will he located at (xq /» Q ). 

J • O } O 





Since there are no perturbation velocities ahead of the Mach cone centered 
on the (x, y, z) system we will define a new variable I = X~Ap. On a 
cylinder with large p the values of 5 will be of order one in the region where the 
Mach cone intersects the cylinder. 
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These functions must be evaluated at the start and end of the vortex strip 
(linearly varying line doublet). 
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To find the interference wave drag between two strips we must integrate the 
product of u from one strip times Vr from the other strip and multiply by the 
element of area over the infinite cylinder. Assume 
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strip 2 begins at 
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with corresponding values for ends 2, 3, and 4 
The wave drag is evaluated by integrating 
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Considering only a single term of the above and letting 
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Therefore, 
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The integrals over 0 may be evaluated as follows 
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The interference wave drag between two strips is 
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the interference wave drag between two elementary horseshoe vorticies is 
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This is the result given on page 35* 
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Vortex Drag.- To evaluate the vortex drag the following integral must 
be evaluated in the Trefftz plane. 
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we can write 




Then, using Stokes theorem, the surface integral may be reduced to a series 
of line integrals around the singularities: 
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The velocities induced by an elementary horseshoe vortex in the 
Trefftz plane are found by taking the limit as x — - «> 
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Ft»r two elementary horseshoe vorticies (doublets) the velocity potential is 
the sum of the individual velocity potentials 
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The interference drag between doublet ) , and doublet 3 is given by 
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Hie contour c will be chosen to be a small circle around each doublet. 
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To evaluate 
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The vortex drag is thus 
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which is the result given on page 35 • 
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RESULTS 


Analysis of the spanwise variation of leading edge thrust using the 
linearly varying vortex finite element code of Appendix E is presented in 
this section. The numerical solution has been compared, to analytic results, 
where available. 

The calculated spanwise distribution of thrust for a highly swept cranked 
wing at incompressible speeds is presented on figure 2. A strong suction 
peak exists at the leading edge discontinuity. Weakening of the singularity 
requires cambering of the wing in this region or elimination of the break in , 
the planform. 

The calculated spanwise distribution of thrust for an aspect ratio 2.5 
wing with a leading edge sweep of 63 . and a taper ratio of .23 is presented 
on figure 3 -fpr several Mach numbers . Transonic and supersonic 

conditions at an angle of attack of 1 radian are shown. The M. =2.0 case 
corresponds to a 3*^° subsonic leading edge. A comparison of the numerical 
result with the exact analytic solution 11 is presented at supersonic speeds. 
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To illustrate the entirely different leading edge suction distribution 
on a planform with forward sweep, the results for the same planform in 
reverse flow (referred to as M M = -0.90) are presented in the same figure. 
Although the suction values obtained by the program are qualitatively 
correct, better agreement is desired. The descrepancies are felt to be due 
to the problems associated with the pressure distributions obtained with the 
spanwise linearly varying panel. In most cases the boundary conditions can 
be satisfied exactly at all control points only at the expense of a considerable 
degree of instability in the spanwise pressure distribution. 

The first attempt to alleviate the spanwise instability in the 
pressure distribution involved the use of spanwise constraint functions. 

This method sacrificed the boundary condition at the outer control point. 

The results, although they were a. vast improvement, were still not 
satisfactory. Thus, additional methods were tried. 
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The first method involves using the pressure distributions obtained 
from a constant vorticity finite element. Since the calculation of the near 
field drag requires only a pressure distribution to be given, no theoretical 
problems arise. The second method uses the linearly varying finite element to 
obtain the pressure distributions, but places the control point at a small 
fraction of the spanwise extent of each panel. The pressures are not strongly 
dependent on the spanwise location and a value of 0.15 was found to give 
satisfactory results. 

There still remain some problems with each method. For planforms with 
large sweeps near the kink, the constant pressure panel gives pressure 
distributions which are felt to be too high near the leading edge. This 
results in negative values of computed leading edge suction in this region. 

The linearly varying panel, with the control point at 0.15 span, yields 
unsatisfactory pressure distributions near cranks. This problem can be 
overcome somewhat by moving the control point nearest the crank to a 
point near the midpoint of the panel. 

At present the program allows the user to select either constant pressure 
panels or linearly varying panels to obtain the pressure distributions. For 
the linearly varying panels the fraction of the span for the control point 
may be specified and this may be overridden at any span station by specifying 
the spanwise location of the control point. Typical results for each method, 
when used on the planform shown in figure 3, are presented as the sample cases 
in Appendix E. 
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CONCLUSION 


The technique presented in this report permits for the first time the 
exact analytical calculation of near field drag at all Mach numbers under the 
constraints of linearized theory. The technique permits the calculation of an 
influence matrix with elements representing the drag which one finite element 
induces on another. The drag may be calculated when the singularity strengths 
are known. The method is applicable to either source or vortex finite elements 
and can be used on planforms of arbitrary shape which may be nonpianar in the 
spanwise direction. 

In addition* the method guarantees the equivalence of the drag as calculated 
in the near and far field, because the form of the influence coefficient 
shows that the near field drag term changes sign as the influenced and 
influencing panel are interchanged. 

Although the technique involves analytic integrations, a large number of 
calculations are required. Since the computer code was not completely 
optimized the computation time required to compute the drag influence matrix 
is on the order of ten times the time required to compute the aerodynamic 
influence matrix. However once calculated the matrix may be stored and need 
not be recalculated if only the twist, camber and thickness distributions 
are changed. 

Additional problems remain with the instability produced in the pressure 
distributions derived with the spanwise linearly varying lifting element. 
Although these problems have been largely overcome by relocating the control 
point or vising pressure distributions obtained from constant pressure panels, 
some problems still remain. 
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APEPDIX J.. 
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APPENDIX B 


TABLE OF INTEGRA IB 
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APPENDIX C 


INTEGRALS FOR THE CALCULATION OF 
MUTUAL INTERFERENCE DRAG 


All of the quantities required to compute the induced drag are composed 
of integrals and derivatives of the basic velocity components induced by 
constant source or constant vortex panels. On pa.ge 11 set of velocities 
was defined such that : 
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Only 5 basic terms occur in the expressions for the velocity components 
for the non-planar source and vortex panels. The functions will be denoted 
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All of the expressions required for the calculation of drag are integrals 
of these functions. If we define 
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then all of the expressions which are required are such that i+j 5 4 


The integrals and derivatives were computed analytically and the results are 
presented in this Appendix. These analytical results were checked by performing 
the differentiations numerically in a computer program. Thus the velocity 
components corresponding to a vortex panel could be- constructed from: 
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Using the derivatives with respect to T, given in Appendix A, we can show 
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These relations were all verified numerically. 
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RELATIONS BETWEEN INTEGRALS 
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APPENDIX D 


SPANWISE INTEGRATION OF NONPIANAR RJNCTIOIB 


It Is desired to integrate expressions of the form 
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Since all of the desired integrals are related to five basic functions, 
the integration can be facilitated with the following scheme. 

Let 

Oo 

,U) 

(!,1 

• o J 

, <S) 

and 



^ O.) 

^ (. 5 1 l B) 

35 Lj J J J 



, a.) , 

f C5,7, S.O) 
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Theh, using the relations given in appendix C 
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However. the following relations exist between the functions. 
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By using the proper choice for a., , a . 3 , a. f 

in each row of the inatricies g. and a different and more useful pair 

of c* and H inatricies may be derived. 
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With the changes we can write 
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A linear trans format ion exists between the vectors 
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However, a trass format ion also exists between A and A and also between h"« , 
H'and h"g,h*' 
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Then- the 


following relations 


may be used to perform the 


required integrations. 
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POTENTIAL FORM DRAG PROGRAM 


PROGRAM INFORMATION 

91 

INPUT DATA 

92 

SAMPLE CASE 

95 


PROGRAM INFORMATION 

A computer code has been written to perform the calculations needed 
to apply the method developed in this report to arbitrary wings. The 
code allows the user a number of options through the choice of inputs. 

For example, the pressure distribution can be calculated using the linearly 
varying vorticity panels or using constant pressure panels. Another option 
is in wing geometry; for simple planforms only leading-and t railing-edge 
sweeps, root chord and number of panels are needed as input, but for 
complicated wings the entire arrays of panel corners may be specified. 

The inputs required for the various options are described in the following 
section. 


91 


•»> Q Kf 


INPUT DATA 


ALL DATA EXCEPT THE TITLE CARD IS INITIALLY INPUT INTO A SINGLE 
LARGE ARRAY. EACH DATA CARD HAS THE FORMAT 112, 5F10. 5. THE FIRST 
NUMBER ON EACH CARD IS AN INTEGER GIVING A LOCATION IN THE INPUT 
ARRAY; THE FOLLOWING NUMBERS ON THE CARD SPECIFY THE VALUES TO BE 
INPUT INTO THAT AND THE FOUR IMMEDIATELY SUCCEDING LOCATIONS. 

THESE NUMBERS ARE IN F FORMAT AND MUST INCLUDE A DECIMAL POINT. 
LOCATIONS LEFT BLANK WILL REMAIN UNCHANGED. 

ALL LOCATIONS ARE INITIALLY SET EQUAL TO 0. 

THE LAST CARD IN EACH CASE MUST HAVE A MINUS SIGN IN COLUMN 1. 

PRECEEDING ALL DATA CARDS IN EACH RUN IS A SINGLE DATA CARD 

SPECIFYING THE ARRAY SIZES USED FOR ALL SUCCEEDING CASES. THE 
CARD CONTAINS TWO VALUES AND HAS THE FORMAT DESCRIBED ABOVE. 

THE LOCATION NUMBER (COLS 1-12) MUST BE 1. 

THERE MUST BE A MINUS SIGN IN COLUMN 1 

THE FIRST VALUE MUST EE AS GREAT AS THE MAXIMUM NUMBER OF 
OF BOUNDARY CONDITIONS IN ALL SUBSEQUENT CASES OF THIS RUN. 

THE SECOND VALUE MUST BE AS GREAT AS THE MAXIMUM NUMBER OF 
PANELS IN ALL SUBSEQUENT CASES. 


THE FIRST CARD IN EACH CASE CONTAINS THE TITLE 


LOC VARIABLE 


DESCRIPTION 


1 


N 
M 

MC 

3 NDRG 

6 SYM 

7 I NFL 


SET THIS LOCATION TO -1. AFTER THE LAST CASE TO 
TERMINATE THE RUN 

THE NUMBER OF PANELS IN THE CHORDWISE DIRECTION 
THE NUMBER OF PANELS IN THE SPANWISE DIRECTION 
THE NUMBER OF CONTROL POINTS IN THE SPANWISE DIRECTION 
IF 0. IT WILL BE SET EQUAL TO M 
=0. THE NEAR FIELD DRAG IS COMPUTED 
<0. NO NEAR FIELD DRAG IS COMPUTED 
0. SYMMETRY ABOUT Y=0. IS ASSUMED 
. NE. 0. NO SYMMETRY 

2 COMPUTE BOTH AERO AND DRAG MATRICIES AND STORE 
ON UNIT 11 

1 COMPUTE NEW AERO MATRIX, READ DRAG MATRIX FROM 
UNIT 11, THEN STORE BOTH ON UNIT 11 
0 COMPUTE NEW AERO AND DRAG MATRICIES 
-1 USE AERO AND DRAG MATRICIES FROM THE PREVIOUS CASE 
-2 READ AERO AND DRAG MATRICIES FROM UNIT 11 
-3 COMPUTE AERO MATRIX AND READ DRAG MATRIX FROM 
UNIT 11 


92 



12 

100 

16 

IWPRNT 

19 

INTMED 

20 

INFPRT 

22 

ND 

23 

MD 

101 

CENT 

102 

SPAN 

103 

ALPHA 

105 

MACH 

108 

CAVG 

109 

SREF 

110 

SWEEPL 

111 

SWEEPT 

112 

ROOT 


113 

RATIO 

114 

RAT I OY 

117 


118 


119 


121 

Z(J) 


IF 0. LINEARLY VARYING PANELS ARE USED TO COMPUTE CP'S 
IF >0 CONSTANT CONSTANT PRESSURE PANELS ARE USED TO 
COMPUTE THE PRESSURE DISTRIBUTION 
IF 0. THE DOWNWASH VALUES ARE NOT PRINTED 
VARIOUS ORDERS OF INTERMEDIATE PRINTOUT (-1. TO 2. ) 

-1. LEAST PRINTOUT 

2. MOST PRINTOUT 

*1 PRINT THE DRAG INFLUENCE MATRIX 
=2 PRINT THE AERO INFLUENCE MATRIX 
THE NUMBER OF CHORDWISE LOCATIONS OF CAMBER INPUT 

THE NUMBER OF SPANWISE LOCATIONS OF CAMBER INPUT 

1. CONTROL POINT AT PANEL CENTROID 

CONTROL POINT AT PANEL CENTER OTHERWISE 
SPAN, ANY CONSISTENT SET OF UNITS MAY BE USED. THIS 
VALUE IS NOT USED IF LOCATIONS 161 TO 161+M ARE USED 

THE ANGLE OF ATTACK IN DEGREES. THIS IS IN ADDITION TO 

ANY CAMBER DISTRIBUTION 
THE MACH NUMBER 

THE REFERENCE CHORD LENGTH ( IF 0. AN AVERAGE CHORD 
THE REFERENCE AREA. IF 0. , SREF=CAVG*SPAN 
THE LEADING EDGE SWEEP IN DEGREES. THIS VALUE IS 
IGNORED IF 112 IS . LE. -1. , WHICH MEANS A PLANFORRM 
SHAPE IS TO BE READ IN. 

THE TRAILING EDGE SWEEP IN DEGREES 

ROOT DIMENSION OR CHORD LENGTH ALONG THE SYMMETRY AXIS 
<0. THE PLANFORM SHAPE IS READ IN FROM 241-320 
=0. THE ROOT DIMENSION IS COMPUTED USING THE SPAN AND 
LEADING AND TRAILING EDGE SWEEPS TO MAKE CAVG=1. 

>0. THIS VALUE IS USED TO CALCULATE THE GEOMETRY 
THE CHORDWISE CONTROL POINT LOCATION. IF 0. DEFAULT IS 
0. 950 FOR MACH > 1. 

0. 875 FOR MACH < 1. 

THE SPANWISE CONTROL POINT LOCATION. IF 0. A VALUE 
BASED ON LOCATION 101 IS USED. 

STARTING CL FOR DRAG POLAR. DEFAULT = 0. 

ENDING CL FOR DRAG POLAR. DEFAULT = 1. 

DELTA CL FOR DRAG POLAR. DEFAULT = . 05 

Z VALUES AT THE Y COORDINATES 
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161 Y( J) Y COORDINATES. (IF ( ( Y(2)-Y( 1 ) )**2+(Z (2)-Z( 1 ) )**2) IS 

0. THE SEMI SPAN IS BROKEN INTO M EQUAL SEGMENTS) 
THERE MUST BE M+l VALUES INPUT 


201 YC ( J) 


241 XL< J) 


281 XT( J) 


Y COORDINATES FOR THE CONTROL POINTS. THESE VALUES 
MUST BE INPUT IF MC > M. IF MC=M, NONZERO VALUES 
WILL BE USED TO OVERRIDE VALUES BASED ON 114 
THE LEADING EDGE COORDINATES AT Y(J) 

THESE VALUES ARE CONSIDERED ONLY IF 112 IS < 0 
XL ( 1 ) CORRESPONDS TO THE COORDINATE ON THE AXIS 
XL (M+l) CORRESPONDS TO THE COORDINATE AT THE TIP 
ANY VALUES WHICH ARE EXACTLY 0. WILL BE CHANGED 
TO MAKE THE EDGE STRAIGHT BETWEEN THE TWO 
SURROUNDING NONZERO VALUES. 

THE TRAILING EDGE COORDINATES. SAME FORMAT AS 241-280 


321 XD( I ) THE X/C VALUES WHERE THE CAMBER IS INPUT. SEE 22 

341 YD( J) THE Y VALUES WHERE THE CAMBER IS INPUT. SEE 23 


401 THE VALUES OF DZ/DX FOR THE CAMBER. 

401 TO 400+ND ARE FOR YD(1) AT XD(1) TO XD(ND) 

401+ND TO 400+ND*2 ARE FOR YD(2) AT XD(1) TO XD(ND) 
THESE VALUES ARE INTERPOLATED TO OBTAIN THE BOUNDARY 
CONDITIONS AT THE CONTROL POINTS. 
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SAMPLE CASES 


This section present samples of the input and output of the computer 
program which is used to calculate the spanwise variation of drag due to lift. 
The planform used for the calculation had a straight leading edge with a 
sweep of 63A35 degrees and a straight trailing edge with a sweep of 45°. 

The Mach number was I.30 and the paneling scheme used 10 panels spanwise 
and 8 panels chordwise. (see figure 3)- The planform has no twist or 
camber and has an angle of attack of one radian. 

There are two cases presented. The first uses spanwise linearly varying 
vorticity panels to compute the pressure distribution on the above planform 
while the second case uses constant pressure panels to compute the pressure 
distribution. The linearly varying panels have the control point placed at 
0.15 of the span while the constant pressure panels have the downwash placed 
at the centroid. The pressures obtained by the constant pressure panels 
are interpolated to the panel edges using a curve fit subroutine. 

A sample input for the first case is presented on pages 97 and 98 - The 
program output begins on page 99- 

Calculated geometric characteristics of the planform are presented on 
pages 100 and 101. The coordinates of the leading and trailing edges 
and the resulting chord dimensions are printed at the y coordinates of the 
panel edges and at the control points. 

Page 102 presents the Cp values at the inside leading edge of each 
panel using linearly varying vorticity panels. 

Page 103 presents the integrated loads at each span station (linearly 
varying panel analysis) as well as integrated loads for the entire planform 
where 

is the span coordinate of the panel midpoint 
is the CL*C/CAVG at the left edge of this row of panels 
is the CL*C/CAVG for this row of panels 

is (x/c) value for the center of pressure at this span station 
is the CD*C/CAVG for this span station 

is the zero suction drag for this span station obtained by 
multiplying the CP value times the control point downwash 

is the leading edge suction for this span station obtained 
from CD0(Y) -CD(Y) 


Y 

CN(Y) 

CL(Y) 

XCP(Y) 

CD(Y) 

CD0(Y) 

CT(Y) 
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Page 104 presents the calculated drag polar. A comparison of the 
calculation and the exact analytic results for the spanwise variation of 
leading edge thrust is presented on figure 3* 

The results of a second sample case using the constant vorticity panel 
analysis on the same planform are presented on pages 105 through 109* 

The downwash presented on page 107 represents the downvash obtained 
using a linearly varying panel influence matrix, with the control point 
located at the centroid, and the pressures obtained using the constant 
pressure panel with the control point also located at the centroid. 
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FORTRAN FIXED 10 DIGIT DECIMAL DATA 



FORTRAN FIXED 10 DIGIT DECIMAL DATA 
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CONTROL POINT AT PANF.L CENTROID 


table of downwash values at control points 
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CONTROL POINT AT PANEL CENTROID 


7500 3.2537 3.1861 0.4370 3.5792 3.1861 -0.39313 
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DRAG POLAR 
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PROGRAM CODE 


The complete program is available from COSMIC (Computer Software 
Management and Information Center, 112 Barrow Hall, University of Georgia, 
Athens, GA 30602) as Potential Form Drag Program, LAR-12236. The program 
subroutines, and their principal functions, are summarized in the following 
table in alphabetical order. 


SUBROUTINE 

AAMAIN 

AIJKL 

AMTX 


FUNCTION 

Controls program flow 

Calculates the vorticity panel influence coefficients 

Reads the main aerodynamic matricies on or off of the 
disk units 


CAMBER 

COUNTV 

DECRD 

DISPLY 

DRGMTX 

DRAG01 

FMTX 

F2F3 

G^MFRY 

HSHLDR 

INTRP3 

LATTCE 

LIFT 

MAIN 

MATRXF 

MPXMLT 

PLNFRM 


Interpolates the input camber to obtain the boundary 
conditions 

Computes the elapsed CPU time for a given program step 
Subroutine used to read the input data 
Prints computed pressures and downwash 

Controls the main flow of the drag influence coefficient 
calculation 

Controls the near field drag influence coefficient 
calculation 

Computes various interpolation coefficients for camber etc. 

Computes functions F2 and F3 in the near field drag influence 
coefficient calculation 

Computes planform geometry 

Linear equation solution routine using a least square 
solution if there are more equations than unknowns 

A third order curve fit subroutine 

Subroutine for planform geometry printout 

Computes planform lift, moments and drag 

Sets up the structure and size of the main arrays 

Routine used to print computed arrays of pressure 
coefficient, velocities, etc. 

Multiplies matricies 

Computes some program geometry 
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SUBROUTINE 


FUNCTION 


POIAR 

RL0GX 

TIMERV. 

TRAPZD 

VTXDRG 

VELCTY 

WVEVTX 


Confutes the 0 and 100 percent suction drag polars 

Computes R and FI between all panel corners for the near 
field drag computation 

Computes the elapsed CPU time for a giver program step 

Obtains the Cp distribution using constrant pressure 
panels 

Computes the vortex drag using a Trefftz plane solution 

Calculates FI, F2 and F3 for the vorticity panel 
influence coefficients 

Calls the routines to perform the drag calculations 
and stores the results on disk units 
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